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The ground state spin of the negatively charged nitrogen-vacancy center in diamond has many 
exciting apphcations in quantum metrology and solid state quantum information processing, includ- 
ing magnetometry, electrometry, quantum memory and quantum optical networks. Each of these 
applications involve the interaction of the spin with some configuration of electric, magnetic and 
strain fields, however, to date there does not exist a detailed model of the spin's interactions with 
such fields, nor an understanding of how the fields infiuence the time-evolution of the spin and its 
relaxation and inhomogeneous dephasing. In this work, a general solution is obtained for the spin in 
any given electric-magnetic-strain field configuration for the first time, and the influence of the fields 
on the evolution of the spin is examined. Thus, this work provides the essential theoretical tools for 
the precise control and modeling of this remarkable spin in its current and future applications. 

PACS numbers: 31.15.xh; 71.70.Ej; 76.30.Mi 



I. INTRODUCTION 

The negatively charged nitrogen- vacancy (NV~) cen- 
ter in diamond has many exceptional properties that are 
highly suited to applications in quantum metrology and 
quantum information processing (QIP). The exciting re- 
cent demonstrations of high precision magnctometryfi"— 
electrometry;^ decohcrence based imaging^"— and spin- 
photoni^ and spin-spini^^— entanglement have each uti- 
lized the ground state spin as a solid state spin qubit. 
Indeed, these demonstrations exploit the interaction of 
the spin with some configuration of electric, magnetic 
and strain fields and the center's remarkable capability 
of optical spin-polarization and readout . The the- 
ory of the ground state spin was presented in part I of 
this paper series^ including the spin's fine and hypcrfine 
structures, and its interactions with the different fields. 
The theory will now be applied to derive the solutions 
and dynamics of the spin for a general field configura- 
tion, including the effects of field inhomogeneities and 
the influence of fields on the spin's relaxation. 

The NV~ center is a point defect of Csu symmetry 
in diamond consisting of a substitutional nitrogen atom 
adjacent to a carbon vacancy (refer to Fig. [ij. The cen- 
ter's electronic structure is summarized in Fig. [2] It con- 
sists of a ground triplet state, an optical excited 
triplet and several dark singlet statesi^ At ambient tem- 
peratures, the fine structures of the ground and excited 
triplet states have single zero-field splittings between the 
TOs = and nis = ±1 spin sub-levels of Dgs ^ 2.87 GHz 
and Des ^ 1-42 GHz respectivelyj^ Zeeman, Stark and 
strain splittings have been observed in the fine structures 
of both triplet states^"— although the Stark and strain 
effects in the ground triplet state are several orders of 




FIG. 1: (color online) Schematic of the nitrogen-vacancy 
center and the adopted coordinate system, depicting: the 
vacancy (transparent), the nearest-neighbor carbon atoms 
to the vacancy (black), the substitutional nitrogen atom 
(brown), the effective magnetic and electric-strain field com- 
ponents (colored arrows), and their corresponding field angles. 



magnitude smaller than in the excited triplet statci^ 

One of the most intriguing properties of the NV~ cen- 
ter is the ability to optically polarize and readout the 
ground state spin^i^l distinguishing the center from 
other paramagnetic defects in diamond, and forming 
the basis for its QIP and quantum metrology applica- 
tions. As depicted in Fig. [2l the process of optical 
spin-polarization occurs due to the presence of a non- 
radiative decay pathway from the excited triplet state 
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FIG. 2: The electronic orbital structure (left) and fine struc- 
ture (right) of the NV^ center at ambient temperatures. The 
observed optical zero phonon line (1.945 eV)2^ and infrared 
zero phonon line (1.190 eV)^ transitions are depicted as solid 
arrows in the orbital structure. The radiative (chain arrows) 
and non-radiative (dashed arrows) pathways that result in 
the optical spin-polarization and readout of the center are 
depicted in the fine structure. Note that the much weaker ra- 
diative and non-radiative transitions that act to reduce spin- 
polarization have not been depicted. 



splittings of the ground triplet state. The microwave 
pulses coherently manipulate the ground state spin and 
result in an optically detectable oscillation in the relative 
population of the = and = ±1 sub-levels. In 
order to optimally control the spin and illicit the maxi- 
mum amount of information and sensitivity from its im- 
plementation in ODMR experiments, a detailed model of 
the time-evolution, relaxation and dcphasing of the spin 
is required. 

In this article, the detailed theoretical construction of 
the ground state spin-Hamiltonian that was conducted 
in Part of this paper series will be applied in or- 
der to produce the solution of the ground state spin and 
its time-evolution in the presence of a general electric- 
magnetic-strain field configuration. The solution will be 
demonstrated by modeling a simple Free Induction De- 
cay (FID) experiment and examining the dependence of 
the FID signal, inhomogcneous dcphasing and spin re- 
laxation on the applied field configuration. This sim- 
ple demonstration will provide insight into the observed 
strong dependence of the spin's inhomogcneous dcphas- 
ing time on the applied fields Furthermore, the spin 
solution will be a useful tool in future applications of 
the spin in quantum metrology and QIP as it clearly de- 
scribes how electric, magnetic and strain fields can be 
used to precisely control this important spin in diamond. 
For example, multimodal decoherence microscopy that 
maps both magnetic and electric noise using the same 
probci^"— 



to the ground triplet state that competes with the op- 
tical decay pathway. The details of the photokinetics of 
the non-radiative pathway are yet to be fully explained, 
but it is believed that the = ±1 sub- levels of the 
excited triplet state are preferentially depopulated and 
the nis = sub-level of the ground triplet state is pref- 
erentially populated, thereby polarizing the population 
into the iris ~ spin state after optical cycling i^i^ The 
preferential non-radiative depopulation of the to^ = ±1 
sub-levels of the excited triplet state also introduces a 
difference in the optical emission intensity between the 
spin sub-levels. This difference in emission intensity can 
be utilized to readout the relative populations of the 
rris = and = ±1 sub- levels of the ground triplet 
state through the measurement of the integrated emis- 
sion intensity upon optical excitationj^i^ 

The center's capability of optical spin-polarization 
and readout enables the implementation of continuous 
wave and pulsed optically detected magnetic resonance 
(ODMR) techniques,20i2i Simple pulsed ODMR tech- 
niques such as free induction decay (FID) and spin 
echo^Sii^ as well as more complicated multi-pulse ODMR 
techniques^'^'^ have been implemented in the center's 
quantum metrology and QIP applications and involve op- 
tical polarization and readout pulses encompassing a se- 
quence of microwave pulses tuned to the fine structure 



II. SOLUTIONS OF THE GROUND STATE SPIN 

The ground state spin-Hamiltonian as derived in Part 
I is 

" ^"^^z + ■ B - -^ExiSl ~ Sy) 

~^J^^yi^xSy + SySx) (f) 

where S — SxX + Syy + SzZ is the total electronic spin 
operator for S = 1, V is the effective spin-spin and axial 
electric-strain field, B is the effective magnetic field, and 
£x and £y are the effective non-axial electric-strain field 
components. Refer to Part I for the explicit definitions of 
the effective fields. As discussed in Part I, the solutions 
of Hgs describe both the interactions of the electronic 
spin-orbit states in the high field limit, where the field 
induced shifts are much larger than the spin's hyperfine 
structure, and also the interactions of the mj = sub- 
set of hyperfine states in the weak field limit, where the 
field induced shifts are comparable to the spin's hyperfine 
structure. 

It is convenient to define the field spin states 
{|0), |— ), 1+)} in terms of the Sz eigenstates {\S,ms)} 
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where tan^g = £y/£x, tan0 = E^/B^ and £±^ = 
£^ + £y. The matrix representation of Hgg in the basis 
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where (f> = 20/3 + (j)£, tan^ij = By/Bx, TZ = + , rections to the energies also become important, such that 



= cos : 



and 



sm ■ 



Therefore, if Bj_ = e'^^^^ = -ATI and 



Bl+B"^ <^ V and B^ /V < TZ, the field spin states 

are approximate eigenstates of Hgs with energies iJg = 
and E± = I? ± 7?.. In this weak non-axial magnetic field 
limit, the spin eigenstates are completely characterized 
by the field angles 0^ and 9 derived from the axial mag- 
netic and non-axial electric-strain field components and 
the energy splitting of the |±) states is governed by the 
magnitudes of the same field components. 

When B\I'D ^ TZ, first-order corrections to the field 
spin states become important, such that 
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where A = B\/2n'D, and second-order perturbation cor- 
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Consequently, in the strong non-axial magnetic field 
limit, the spin eigenstates and energies become depen- 
dent on the non-axial magnetic field direction 05 and 
the dimcnsionless ratio of the field magnitudes A. 

Figurc[3]contains polar plots of the dimcnsionless split- 

ting parameter (l — 2A sin cose/) -I- A^) ^ of the |±) spin 
states as a function of the azimuthal 4) and polar field 
angles. Clearly the most interesting field configurations 
have the parameters A ~ 1 and 9 ^ since for these 
parameters the energies depend sensitively on the field 
angles. Such a field configuration was used in the recent 
electric field sensing demonstratior^ to sensitively detect 
both the magnitude and orientation of an applied electric 
field. Future implementations of the ground state spin as 
a field sensor should seek to exploit these high sensitivity 
field configurations. 

The matrix representation of the interaction of the 
ground state spin with an oscillating microwave magnetic 
field M in the basis {|0), |-), |+)} is 
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where Ai± = 



Ml+Ml, (pm = 2(j)M + (t)£ and tan = A^^/TW a;- The off-diagonals of the above ma- 
trix representation indicate that the transitions between 



(a) (b) 

FIG. 3: (color online) Plots of the dimensionless energy splitting parameter (l — 2A sin cos + A^) ^ of the |±) spin states as: 
(a) a function of the azimuthal < < 27r and polar < 6^ < tt field angles for A = 1; and (b) as a function of for — ^ and 
different values of A as indicated. Coordinate axes are provided for reference to figure [1] where the field angles — 2(f>B + 4'e, 
4>e, 4'B and 9 are defined. 



the field spin states induced by the oscillating microwave 
field are also dependent on the static field angles. 

Using Fermi's golden rule,— the transition rates Wo^± 
between the |0) and |±) spin states correct to zero-order 
in the static non-axial magnetic field B± are proportional 
to the absolute square of the off-diagonal elements, such 
that 

Wo^± (X Tsin6lcos(^™) (6) 

As the transition rates to the |±) spin states depend dif- 
ferently on the static fields for a given microwave po- 
larization, the transitions to the |±) spin states can be 
controlled via the static fields, or conversely for a given 
static field configuration, the transitions can be con- 
trolled via the microwave polarization. Figure 0] depicts 
the dependence of the transition rates on the microwave 
polarization and static field configuration. As can be 
seen, the individual transitions can be selectively ex- 
cited using orthogonal non-axial microwave polarizations 
ipM — + f when 6* ^ f • Hence, linearly polarized 

microwaves in conjunction with static field control or (as 
has been previously demonstrated)^ circularly polarized 
microwaves can be used to selectively excite individual 
spin transition in situations where the splitting of the 
|±) spin states (= 27?.) is too small to selectively excite 
the individual transitions using microwave frequency se- 
lection alone. 
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FIG. 4: (color online) Plots of the normalized transition 
rates VKo->- oc ^(1 -|- sin 6* cos 0m ) (blue) and Wo-^+ oc 
i(f — sin 6* cos 0m) (red) between the |0) and |— ) and the |0) 
and |-|-) spin states, respectively, as functions of the non-axial 
microwave polarization (j>M , given (j>e = and different values 
of the static field angle 9. Note that the transition rates are 
equal for all values of (pM when ^ = 0. 



III. TIME-EVOLUTION OF THE SPIN: 
SIMULATION 



FID 



Assuming the weak static non-axial magnetic field 
limit (i3_L < V and B\lV < 7^), which typically occurs 
in most of the current applications of the ground state 
spin, the spin solutions can be used to accurately model 
the time-evolution of the spin in a simple FID ODMR 
experiment. As depicted in Fig. [Sj the FID sequence 
is comprised of optical pulses that polarize the spin at 
t ~ ~tr and readout the spin at i = i^, as well as mi- 
crowave 7r/2 pulses that coherently manipulate the spin 
before and after the period of free evolution r. Note that 
the MW pulse sequence depicted in Fig. [S] differs from 
the conventional ESR sequence by the final 7r/2 pulse, 
which projects the accumulated phase into a population 
difference between the nis = and rus = ±1 spin sub- 
levels. 

In the model of the FID experiment, the static fields 
will be considered to differ during the period of free evo- 
lution, such that before and after the period of free evo- 
lution the static fields are described by the parameters 
{V, TZ, 9, 0f ) and the |±) spin states have energies haj±, 
whilst during the period of free evolution the static fields 
are described by the parameters {V , TZ' , 9', (j)'^) and 
the |±) spin states have energies huj'j.. For simplicity, the 
changes in the static field configurations are assumed to 
be adiabatic and infinitely sharp at t = ±r/2 and the 
microwave field is assumed to selectively excite the tran- 
sitions between the |0) and |— ) spin states with a tuned 
microwave frequency w ~ uj- . 

This model FID experiment is a generalization of the 
FID experiments that were conducted in the spin's re- 
cent electric field sensing demonstration^ and one of the 
spin's magnetic field sensing demonstrations^ and will 
consequently describe how the combined effects of elec- 
tric, magnetic and strain fields will influence such sensing 
demonstrations. However, note that the objective of pre- 
senting this model FID experiment in this work is not to 
discuss sensing techniques, but to provide the necessary 
theoretical details to discuss the effects of inhomogeneous 
fields and lattice interactions on the relaxation and de- 
phasing of the spin in the following sections. 

The state of the spin at a given time during the FID se- 
quence \t) can be written in terms of the spin eigenstates 
of the static fields at that time 



c^{t)\-)+c+{t)\ + 

^cJ_{t)\-)' + c'4t)y 
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where the coefficients Ci{t) and c^(t) (i = 0,— ,-f) arc 
related at i = ±t/2 by the basis transformation T : 
{|0), 1+)} ^ {|0)', h)', 1+)'} given by the matrix 
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FIG. 5: (color online) The free induction decay (FID) se- 
quence of optical pulses that polarize the spin at t = —tr 
and readout the spin at t = tr, as well as microwave ir/2 
pulses that coherently manipulate the spin before and af- 
ter the period of free evolution r. The static fields and 
spin state energies differ during the period of free evolution 
{T>' , TZ' , 6' , (p'g, &jj-) compared to before and after the pe- 
riod of free evolution {T>, TZ, 6, (fis, fkjj±). The microwave 
field is assumed to selectively excite the transitions between 
the |0) and |— ) spin states with frequency lu. 
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(8) 

where S4>e = ^{4>'e - M and d9± = ^{9' ±9). 

Introducing the density operator p(t) = \t){t\, the ma- 
trix representation of p{t) in the spin eigenstate basis of 
the fields at t ({|0), |-), |+)} or {|0)', |-)', |+)'}) is 



rp(i) -tr<t<-^. 



<t <U 
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where the elements of the matrices p{t) and p'{t) are 
Pij{t) = c,{t)c*{t) and p',^{t) = c'i{t)c'*{t) respectively, 
and at i = ±-| the matrices are related by the transfor- 
mations p'(-§) = r-^p{-^)r and p(§) = Tp'{^)r-\ 

In order to simplify the treatment of the interaction 
of the spin with the microwave pulses, the rotating wave 
approximation can be adopted and the density operator 
transformed into the rotating reference frame, such that 
pit) ^ p{t) and p'{t) p'it), where p^^it) = po±{t) = 

Po±{t)e-^"\/Ut) = Po±W = Po±(i)e-'"*, and p,,(i) = 
Pij{t) and p'jj{t) ~ p'ij{t) for the other density matrix 
elements.— The effect of the 7r/2 microwave pulses can 

and 



therefore be described by p{- 
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p{tr) 

trix representation 



'■^2 p(^)e*'^ 2 ^ where using the microwave ma- 
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(10) 



and (JtiL is 



the duration of the tt/2 pulse. 

Due to the process of optical spin-polarization, the first 
optical pulse of the FID sequence will incoherently polar- 
ize the ground state spin such that at i = —tr, 



Pn 
p{~U) - I p_ 
p+ 



(11) 



where pi are the populations of the spin sub-levels and 
Po +P- +P+ = 1- Note that since it is believed that the 
spin-polarization process does not discriminate between 
the population of the vis = ±1 sub-levels, it is expected 
that p+ p_. The effect of the first tt/2 microwave pulse 
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where Sp = pq ~ p_ and the simplifying assumption 
P-i- ~ P- has been made. After the sudden change in 
field configuration at t = — ^. p(— ^) is transformed into 
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where = cos (5(/)£ cos — i sin cos (5(?-|- and 

7^_i_ = — cos sin (50_ -\- i sin 6(j)£ sin 69^ are elements 
of the basis transformation matrix. 

The coherent time evolution of the ground state spin 
during the free evolution period is governed by the 

Schrodinger equation -^p{t) = p{t),H'gg ^ However, 

the spin also interacts with incoherent time-dependent 
perturbations such as crystal vibrations, the thermal ra- 
diation field and fluctuating fields from local magnetic 
and electric impurities during the free evolution period. 
These incoherent perturbations induce transitions be- 
tween the spin states that lead to spin relaxation "fl^ and 
dephasing rates jfp which are characterized by the Ti 
and T2 times of the spin, respectively)^ Accounting for 
both the coherent and incoherent evolution of the spin,— 
the following matrix equation is obtained 



+Y-oP~- + 7+0P++ 
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where 5uj'+ 



Lo and Slo'j^, 



Note that 



the coherent coupling between the NV spin and proximal 
impurity spins will not be considered in this work. 

The solutions of the off-diagonal elements at i = ^ are 
simply 

P'UI) = PWi\) = P;,+ (-^)e-^+o-e'^"+- 



p:_(-) = p'_+(-)=p'_+(--)e 



(15) 



whilst the solutions of the diagonal elements p'ooi^), 
/5'__(§) and p^^_(§) are more complicated, but can still 
be obtained analytically as combinations of exponential 
functions of the relaxation rates. It follows that after 
the second sudden change in the static fields at i = ^, 
p{^) = T p' {^)T~^ , and after the second 7r/2 microwave 



pulse at t = tr, p{tr) = e~*'^5p(^)e''^5 . The second 
optical pulse at t = tr reads out the proportion of the 
population in the TOs = sub-level, such that the optical 
emission intensity /(r) oc pQo{tr), where 

I [poo(^) + |71_|V-(^) + \T-+\'pWil) 
-y [\T—\^e-''-«-' cos 6lj'_t 
+ \T-+\'^ cos 6uj'_^T 
-\r—\^\r-+\^e'^+-'cos6uj'_^_T 



POo{tr) 



-Nir) 



6p 



0{t) 
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The first term N{t) is not oscillatory and depends 

only on the diagonal elements Poo(~i)i P' ^'^'l 

p(|_^(— ^) and the relaxation rates 7^^. The second term 
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0(r) has oscillatory components with frequencies cor- 
responding to the different frequency shifts and 
Suj'_^ ) and the contributions of each oscillatory compo- 
nent are dependent on the state couplings (|7t | and 

|71_l-|) and the dephasing rates 7^ . The oscillatory term 
therefore offers a great deal of information about the spin 
eigenstates and their energies during the period of free 
evolution. Furthermore, since the observed Ti and T2 
times of the ground state spin typically differ by at least 
on order of magnitude the change in the non-oscillatory 
term over the lifetime of the oscillatory term is negligible, 
and thus can be effectively ignored in an observation of 
the oscillatory term. 

The oscillatory term can be observed by conducting 
ODMR measurements of an ensemble of spins or con- 
ducting many ODMR measurements of a single spin. For 
a given measurement of an ensemble of spins, the field 
parameters (T>,TZ,d,(j)£) and {V ,TZ' ,9' , (j)'^-) will poten- 
tially differ for each spin within the ensemble due to in- 
homogeneities in the fields. Likewise, for an ensemble 
of measurements of a single spin, the field parameters 
will potentially differ between each measurement due to 
differences in the preparation of the spin and the fields. 
These ensemble inhomogeneities introduce an additional 
dephasing decay in the observation of the oscillatory term 
and can be accounted for by introducing statistical distri- 
bution functions of the field parameters and calculating 
the expectation value {0{t))P 



IV. THE EFFECTS OF INHOMOGENEOUS 
FIELDS 

The total dephasing rate of the spin due to interac- 
tions with incoherent time-dependent fields and inhomo- 
geneous static fields is characterized by the Tj* time of the 
spin<^ In the recent electric field sensing demonstration;^ 
it was observed that the Tj* time of the ground state 
spin was highly dependent on the field configuration, 
such that it obtained a maximum in the absence of an 
axial magnetic field and sharply decreased as the axial 
magnetic field was increased at a rate that was inversely 
related to the non-axial electric-strain field. This ob- 
servation highlighted the significant influence that the 
static fields have on the dephasing of the spin and the 
potential to control the susceptability of the spin to dif- 
ferent noise sources. The dephasing due to static field 
inhomogenieites will be discussed in this section and the 
dependence of the incoherent dephasing rates 7^ on the 
static fields will be discussed in the next section. 

Since the ground state spin interacts very weakly with 
axial electric-strain fields, the effect of the variation of T) 
and T)' between the measurements of a single spin will be 
negligible compared to the variations in the other field 
parameters. Note that this is not necessarily the case 
for an ensemble of spins because local strain fields can 
vary significantly between lattice sites. Nor is it nec- 
essarily the case for single spins or ensembles of spins 



if temperature varies appreciably during the conduct of 
the measurements 1^ Nevertheless, the variations in V 
and V are considered negligible in the following. It is 
reasonable to expect that the field parameters TZ and 
TZ' will have statistical distributions f{TZ/h, /i-^,, a-jz) and 
f{TZ'/h,fiTi>,a-jzi) that have the same distribution func- 
tion /, but different mean values fi-ji and fi-jii and differ- 
ent variances a-ji and a-jz' in units of frequency. 

The distribution functions of the frequency shifts Slu'^. 

and Suj'^ can be constructed using the distributions of 

7^ and W asia 



Jo 




^+-('5w+_,Ai+-,cr+_) = -f{-5uj'^_,iin',<yn') (17) 



where the explicit expressions for the means /i± and /i-j 

and the variances a± and a-\ in terms of /z-r., ct-r. 

and a-ji' depend on the distribution /. For example, if / is 
the normal distribution, then the expressions are simply 

and = 2aTi,' ^ 

Using the distributions of the frequency shifts, the ex- 
pectation value of the FID oscillatory term for an ensem- 
ble of measurements of a single spin is 



(0(r)) 

/OO 
cos 5uj'^ tF_ {Slu'_ , , )dSuj'_ 

/OO 
cos diu'_^TF+ (<5w; , pi+ , a+)d6w'+ 
'OO 

/OO 1 
^-cos<5c.;_r 

F+.{^5uj'+_,fi+^,a+-)dSu'^_ (18) 



Note that the expectation values of the state couplings 
involve just the field angles d(p£ and 69±. 

The above expression demonstrates that (0(t)) is po- 
tentially complicated for general state couplings and dis- 
tribution functions and that it is difficult to extract all 
of the information encoded in the oscillatory term. A 
clearer analysis of the oscillatory term is obtained by per- 
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forming a Fourier cosine transformation, 



Af{Bz/h, ^B,(7i3) by a normal distribution Af(x,fi,(7) 



(OH) 

2 



(0(r)) cos lyrdr 



J —OO 

/oo 
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where L(iy,x,'-f) = - f , : P — vi H — 2 , P , — vr I is a sum 
of Lorentzian distributions. Consequently, it is clear 
that {0{v)) is comprised of a collection of lines at 

V =±yu_, ±^4_, ±/i_| which have composite lineshapes 

of L and the distribution functions of frequency shifts. 
The Fourier analysis of the oscillatory term therefore pro- 
vides the frequency shifts, the state couplings as well as 
information about the distributions of the field parame- 
ters from the locations, intensities and shapes of the lines. 
This additional information encoded in the lineshapes 
can be used to infer details about the statistics of the 
local environment of the spin, a notion which (through a 
different approach) forms the basis of the recent propos- 
als of decohcrence imaging)^"— 

The distribution function / of the field parameters TZ 
and TV can be itself constructed from the distributions of 
the electric-strain £± and magnetic field components. 
Let the distributions of the electric-strain and magnetic 
field components be e and /3 respectively, then^ 

-T^ / , e(^,Ai£,cr£)p(^,MS,cre)a-r-d— 

dTl Jn>^Bl+£l n, h h h 



d f , .^/Bz I £± Bz 

-T^ / , <^[^,f^£,<^£)P[^,tJ-B,<^B)d—d— 

dW Jn'>y/Bfh£i ^ ^ h h 

(20) 

where it has been assumed that the variances of the 
electric-strain and magnetic field components are inde- 
pendent of their mean values. 

The above construction can be demonstrated using the 
simple case where the variance of the non- axial electric- 
strain field ug is negligible compared to the variance 
of the axial magnetic field and the mean values of 
the field components ^£ and /ig. Due to the domi- 
nance of paramagnetic impurities over electric impuri- 
ties in diamond,— this simple case is applicable to most 
applications of the ground state spin. Modeling the 
electric-strain field distribution e = 5(£±/h — jJ-e) by a 
delta function and the magnetic field distribution /3 = 



e 2<T^ /v^TTCT^, the distribution function of TZ is 







u <1 



(21) 



u> 1 



where u = Tl/h^e, fi = Hs/f^E, = crgZ/if: and N-jzifi, cr) 
is a normalization constant. Note that an analogous ex- 
pression can be obtained for TZ' by substituting the re- 
spective mean values and variances of the field compo- 
nents. 

Figure m (a) contains plots of the above distribution of 
TZ for the relative magnetic field variance <t = 0.1 and 
for different relative axial magnetic field mean values /i. 
As demonstrated by figure [5] (b) and (c), the mean of 
the distribution of 7?. varies as nn/ i^-s = y^^M-l, which 
would be expected from the relationship TZ ~ \/B1 + S\, 
and that the variance ctti/ ^is depends sensitively on the 
relative axial magnetic field mean ^, except for /i 3> 1, 
where the variance becomes approximately independent 
of /i. 

Consider a FID experiment where the static fields do 
not differ in the period of free evolution. Such FID ex- 
periments are typically used to measure the T2 time of 
the ground state spin via the width of the single line 
that occurs at 6uj'_ = in the Fourier spectrum of the 
oscillatory term. Noting that for such an experiment 
{\T-+?) = (IT— pr.+ P) = 0, = <Jn = on' and 
8T) ~ 0, the expression for the Fourier spectrum 
simplifies to 

/OQ 
L{v, Slj'_,Y-o)F- (S^'- ' , cr^)dSuj'_ (22) 
'OO 

where 

F_(<5t^^,Ai-,(7-) = 

f{TZ/h, fin, on)f{5i^'- + 1Z/h, fin, (Jn)dTZ/h 

(23) 

and assuming the simple case of negligible electric-strain 
field variance, the distribution / is given by (HI]). Since 
the lineshape is a composition of L and F- , the width of 
the line, and thus the T2* time of the spin, will depend on 
both the dephasing rate 7^q and the electric-strain and 
magnetic field distribution parameters that form F-. 

The distribution F_ is plotted in Fig. [7] for differ- 
ent values of the relative axial magnetic field mean /i 
and variance a. The plots clearly demonstrate that the 
width of the F- distribution is highly dependent on fi, 
such that it increases significantly for small increases in 
/i until ^ ~ 1. In the limit ^ 3> 1 it can be seen that 
the distribution has reached its maximum width and is 
very similar to a normal distribution. This is consis- 
tent with the distribution of TZ being dominated by the 
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(a) 




0.000 



FIG. 6: (color online) The distribution flJZ/ h, ^tz, (Jn) con- 
structed for the simple case where the distribution of the non- 
axial electric-strain component has negligible variance com- 
pared to the variance of the axial magnetic field component 
(jg and the mean values of the field components and /^g. 
(a) plots of the distribution as a function of w = TZ/hfie for 
fJ- = I^B/fJ-e = 0, 1/2, 1, 3/2, 2 (in sequential order left to right) 
and the same variance a = ctb/ fJ-e = 0.1. (b) and (c) are plots 
of the mean p,Ti/fJ-e and variance a-jz/fie of the distribution as 
a function of the axial magnetic field mean fi given a = 0.1. 



distribution of Bz when the mean axial magnetic field 
is much larger than the mean non-axial electric-strain 
field. Hence, it can be concluded that due to the dom- 
inance of magnetic inhomogeneities, the contribution to 
the spin's Tj* time from statistical inhomogeneities in the 
static fields is dramatically reduced for field configura- 
tions where ^ — fJ^ef l^s < 1- This conclusion is consistent 



with the observations of the recent electric field sensing 
demonstration^ and has significant implications for the 
spin's other sensing and QIP applications. 
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FIG. 7: (color online) Plots of the distribution 
F-{5(jjL, fi-,a-) that describes the contribution of static 
field inhomogeneities to the shape of the line that occurs 
in the FID experiment that is used to measure the T2 time 
of the ground state spin, (a) the distribution for different 
values of /X = fiB/fJ-e ~ 0,1/2,1,5/4 and a — as/f^e ~ 0.1. 
(b) the distribution for different values of a = 0.1, 1/2, 1, 3/2 
and ^ = 0. Note that each distribution has been normalized 
such that its maximum is 1. 



V. SPIN RELAXATION AND DEPHASING 

As noted in section III, the ground state spin interacts 
with time-dependent incoherent electric and magnetic 
fields and crystal vibrations, which introduce the spin re- 
laxation "fl^ and dephasing 7,^ rates into the evolution of 
the spin. The incoherent electric and magnetic fields arise 
from the thermal radiation field and fluctuating magnetic 
and electric impurities and the crystal vibrations arise 
from the thermal motion of the crystal nuclei. Typical of 
spin systems, spontaneous radiative emission and contri- 
butions from the thermal field are negligible and can be 
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safely ignored)^ As identified in Part I, the spin's inter- 
action with electric fields is much smaller than its interac- 
tions with magnetic and strain fields. Consequently, the 
relaxation and dephasing rates are expected to be well 
described by just the spin's interactions with magnetic 
impurities and lattice vibrations. The contributions to 
relaxation and dephasing arising from interactions with 
magnetic impurities, including their dependence on the 
static magnetic field, has been studied in some detail for 
NV spins in both type lb and type Ila diamond4^"— 
The additional influence of electric-strain fields on these 
interactions has not yet been studied and this will be 
the subject of future work. In this section, the contribu- 
tions to relaxation and dephasing arising from the spin's 
interaction with lattice vibrations will be theoretically 
developed for the first time. 

The linear interaction of the spin with the vibrations 
of the crystal is described by the potential^ 



where as for the interaction of the spin with a static 
strain field?^ Qu,E,q is defined as the u*'' mass-weighted 
normal displacement coordinate of symmetry (E, q) that 
corresponds to an eigenmode of the ground triplet with 
frequency uJu,e in the harmonic approximation, Rq are 
the ground state equilibrium coordinates, and bu,E,q and 
Hj. E q vibration annihilation and creation opera- 

tors respectively. Constructed in an analogous fashion to 
the matrix representation of the spin's interaction with 
the strain field 14tr,~ the matrix representation of Vvih 



in the spin basis {|0), 



-)}is 



vib 



-EE 



u,q 



dVNein,R) 
dQu,E,q 



2uJu,E 



ibu,E,q+bi,E.n) 



(24) 



; <l>e 



.- <l>s 



KM 



' 2 ce 



e 2 ce) 



V2 



( 



2 Co 



e' 2 se — e ' 2 ce) ^(e' 2" ce -I- e 

2 2 2 

— SgCtf,^ —CBC^f: 

s|) -ceCff,^ seCff,^ 



^ Cb) -%(e*2^Ce— e ^ Se)\ 



e 2 se) 



Qu, 



E.,y 



(25) 



where g„ = -V2s2,5{ai\\dVNe/dQu,E\jiJ\e), x = 
_\/^2^+£i^^ s„ ^ are the spin coupling coefficients, 

2V2S2,5 ' 

(aill ||e) are the molecular orbital reduced matrix ele- 
ments of the center (refer to Part I for further details), 

and Qu.E^q = K^E^q + bi.E,q- 

Applying time-dependent perturbation theory, the 
above matrix representation can be used to derive the 
first order spin-lattice transition rates, which will con- 
tribute to both the relaxation and dephasing rates of the 
spin. The first order spin-lattice transition rates are 



W. 



±^0 



W, 



vib(i) _ T^x q'^{^±)pE{uj±) 
h uj± 

TTX^ q'^{uj±)pE{uJ±) 



vib(l) 
O-S-ib 



h 



[nTioJ±) + 1] 
nT{uj±) 



h wj 



Wl':^i' = - cos' 6 
h 



UJ+- 



nTiuj+^) (26) 



where q'^{uj) is the average of over all E symmet- 
ric vibrations of frequency w, pe{^) is the density of 
E symmetric vibrations at frequency w, and nx{i-tj) = 
^^^^hu)/kBT _ 2) is the mean occupation number of ther- 
mal vibrations given by the Bose-Einstein factor. 

The NV transition frequencies a;± and occupy 

the very low frequency end of the vibrational frequen- 
cies of diamond, which range from zero to the Debye 
frequency of diamond ud = 38.76 THz4^ In the low fre- 
quency limit, the long wavelength acoustic modes of the 
lattice have the well known Debye density Pe{'^) ~ Pa^"^ ^ 
where pa is a constant related to the acoustic velocity 
in diamond, and the electron-vibration interaction for 
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the non-local acoustic modes has the approximate form 
q'^{uj) wg^, where is a constant.— Additionally, for 
temperatures UbT ^ fy^±, fv^^ , the thermal occupa- 
tions of the vibrational modes with frequencies corre- 
sponding to the NV transition frequencies can be approx- 
imated by nxiuj) ~ ksT/huiM- The first order transition 
rates for temperatures fcsT 3> huj±, fno^ are thus ap- 
proximately linear in T 



vib(l) 
±->-0 
b(l) 



w, 



vib(l) 



b(l) 



Aa cof? 0ujI_T (27) 



where Aa = T^ksPaQa/^'^- Noting that for an axially 
aligned magnetic field, cos^ 9 = ABl and uj± = 
{V zLTZ)^, it is clear that the first order transitions be- 
tween different spin states depend differently on the field 
parameters. 

Since the electron-vibration interaction q^{ui) and the 
density of vibrational modes psi^jj) increase at higher 
vibrational frequencies, the first order transitions will 
only be the dominant spin-lattice mechanisms at low 
temperatures, where there is only appreciable occupa- 
tions of the low frequency vibrational modes. At higher 
temperatures, the occupation of the more numerous and 
strongly interacting higher frequency modes ensures that 
elastic and inelastic Raman scattering of vibrations will 
become the dominant spin-lattice mechanisms4i The in- 
elastic scatterings will contribute to both relaxation and 
dephasing, whereas the elastic scatterings will contribute 
to just dephasing. Note that as the two-vibration absorp- 
tion/ emission transitions involve vibrations of even lower 
frequencies than the first order transitions, they will be 
negligible at all temperatures. Expecting the most sig- 
nificant contributions to be from vibrational modes with 
frequencies uj 3> a;±, w-i , the elastic and inelastic Ra- 
man scattering rates are approximately 



W. 



ib(2) 



±->0 



ib(2) 



0-nt 
■2 



W. 



■ib(2) 



1, 



-Wr 



vib(2) 



0->0 



g2 {uj)p%{u) 



nT{io)[nT{Lo) + V\duj 



W. 



vib(2) 



±-)-± 



^[l + x'(l±sin0cos30£)+x1 



nT{uj)[nT{uj) -f Ijdw 



(28) 



The dependence of the elastic scattering rates of the 
|±) spin states on the static field angles is particularly in- 
teresting. Figure m contains polar plots of the dimension- 
less scattering parameter [1 -I- x^{l ± sin 6* cos Si/)^) -I- x^] 
as a function of the 9 and field angles for different 
values of X- As can be seen, for |x| > the elastic scat- 
tering rates are minimum at (f)£ — 0, mimicking 
the structural symmetry of the defect center. Hence, 
it appears possible for the orientation of the non-axial 



electric-strain field to be tuned so that the elastic scatter- 
ing rates are minimized or maximized with the difference 
in the minimum and maximum rates 2%^ determined by 
the ratio of spin coupling coefficients x- Note that x is 
currently unknown. 




FIG. 8: (color online) Plots of the dimensionless elastic vi- 
bration scattering parameter [1 + X^(l ~ sin cos 3(/>£ ) -f x^] 
as: (a) a function of the azimuthal <</)£< 27r and polar 
Q < < n angles for x — Ij (b) as a function of c^g for 
Q = and different values of x as indicated. Coordinate axes 
are provided for reference to figure [T] where the field angles 
(j)£ and 6 are defined. 



Considering temperatures not so high that optical 
modes are appreciably occupied in thermal equilibrium. 
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there will be two distinct contributions to the integrals 
in the Raman scattering rates. The first will be from the 
acoustic modes which have electron- vibration interaction 
q^{ijj) ~ ioq^ and mode density pe{^) ~ Pa^^^ ■ The sec- 
ond will be from the strongly interacting local modes of 
the NV center which have frequencies u>i ^ 65 meVJ^ 
The contribution from the local modes can be repre- 
sented by the electron-vibration interaction q^icoi) « qf 
and a sharp peak in the density of modes centered at 
Pe{^jJi) — Pi with width ui. Given the separable contri- 
butions, the integrals can be evaluated and the Raman 
scattering rates reduce to 



w: 



vib(2) 



±-»-0 
2 



ib(2) 
-i>± 



w. 



vib(2) 



X 



A^inT{uJi)[nT{ui) + l]+Al 



W. 



ib(2) 



i [1 + X^(1 ±sin6lcos3(/)£) + x''] 
A^inT{ui)[nT{uJi) + 



(29) 



where Af = nqf pfaf/fi^ujf. Note that the integral over 
the acoustic modes was evaluated in the limit — > oo 
in order to obtain the simple factor. Given the tem- 
peratures being considered, for which high frequency op- 
tical modes are not appreciably occupied, this extension 
of the integral is expected to be inconsequential. Hence, 
the Raman scattering rates depend on temperature in 
two distinct ways due to the distinct contributions of a 
few strongly interacting local modes and many weakly 
interacting non-local acoustic modes. 

Combining the magnetic (W^) and spin- lattice 
^^vib(i) ^^^^ ^vib(2)^ contributions, the relaxation and 
dephasing rates of the ground state spin are 



r 

7±o 








_^vib(2) 






±^ 








= 






^^vib(2) 



V 

7±o 



P 

7+- 



^(7io- 



-7S±) + M^o% + W^±-.±+W^( 



b(2) 
0^0 



vib(2) 



ib(2) 



(30) 



For a simple ODMR experiment using the aj_ transi- 
tion, the spin relaxation Ti and dephasing T2 times are 
defined by 1/Ti = jLq + 7o'_ and I/T2 = j^q, which 
using the above are explicitly 



1 
1 



2Tbi + 2TyibiUj'tT + 2rvib2«T(wi)[nT(w;) + 1] 

vibSJ 



1 

^ + rs2 + 



^7 + (3 - sin 61 cos 30f:) + 
X 



[rvib2"-T(t^;)["T(w;) + 1] + rvibs?"^ 



vib3 



where Tvibi = X^A, rvib2 = X^-^fi and 
47rx^k^A^/15h^ are constants that are independent of 
the static fields and temperature, and Tbi and Tb2 are 
the magnetic contributions that are dependent on the 
static fields, but effectively temperature independent. An 
analogous expression can be simply derived for the uj+ 
transition. 

Noting that 71t(w/) [71^(0;;) + 1] « nxii-^i) for ksT < 
hjji , the contribution to 1 /Ti from inelastic Raman scat- 
terings of strongly interacting local modes has been ex- 
perimentally observed 1^ Likewise, the T and con- 
tributions from the weakly interacting acoustic modes 
have also been observed4^ The spin-lattice contribution 
to 1/Ti only depends on the static fields through the 
presence of uPi in the linear temperature term. This 
dependence on the static fields has not yet been ob- 
served, which is most likely due to the insignificance of 
the linear term at ambient temperatures and the fact that 
most previous measurements have been performed us- 
ing NV ensembles, where resonant interactions between 
NV sub-ensembles and between NV centers and PI cen- 
ters at magnetic fields around B ^ Q, 0.96, 1.44 1.68 
GHzf^i^ will most likely have masked the relatively 
weak dependence of the linear term. The dephasing 
rate I/T2 is dominated by the contributions from mag- 
netic interactions with little observed temperature de- 
pendence in small static ficldsi^ Since the magnetic con- 
tribution is governed by the impurity concentration)^ 
it may be possible to observe the spin-lattice contribu- 
tion in highly pure samples. It would indeed be inter- 
esting to observe the tuning of I/T2 using an electric- 
strain field via the spin-lattice elastic scattering parame- 
ter [1 /x^ + (3 - sin 6* cos 30£ ) + x^] ■ 



VI. CONCLUSION 



As described elsewhere, the magnetic contributions are 
highly dependent on the static magnetic field j^i^i^i^i^ 
but are essentially temperature independent for temper- 
atures > 20 K, due to the impurity spins easily reach- 
ing equal Boltzmann populations of their spin sub-levels 
at low temperatures^^ The spin-lattice contributions are 
instead weakly dependent on the static fields, but have 
distinct functions of temperature that arise from different 
interactions with lattice vibrations. 



In this article, a general solution was obtained for the 
NV spin in any given electric-magnetic-strain field con- 
figuration for the first time, and the infiuence of the fields 
on the evolution of the spin was examined. In particular, 
the control of the spin's susceptibility to inhomogeneities 
in the static fields and crystal vibrations was examined 
in detail. The analysis of the effects of inhomogeneous 
fields revealed the field configurations required to switch 
between different noise dominated regimes. The analysis 
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of the spin's interactions with crystal vibrations yielded 
observable effects that are consistent with previous ob- 
servations and also the basis for future investigations into 
the potential tuning of the spin's dephasing rate. Hence, 
this work has provided essential theoretical tools for the 
precise control and modeling of this remarkable spin in 
its current and future quantum metrology and QIP ap- 
plications. 
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